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<]^ Abstract 

The amplitude for the collision of two hadrons with the lowest order pomeron loop is calcu- 
lated. Numerical calculations show that the loop contribution to the amplitude begins to dom- 
inate the single pomeron exchange at rapidities 8 — 10. Full dependence of the triple-pomeron 

rZ^, \ vertex on intermediate conformal weights is taken into account. 
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In the framework of QCD with a large number of colours iV c 3> 1, strong interactions are me- 
diated by the exchange of interacting BFKL pomerons, which split and fuse by triple pomeron 
vertices. This picture can be conveniently described by an effective non-local quantum field 
theory [1]. From the BFKL approach this theory inherits the fundamental property of con- 
formal invariance, which is only broken by interactions with colliding hadrons. In terms of 
the relevant Feynman diagrams contributions standardly separate into tree diagrams and di- 
agrams with pomeron loops, which corresponds to division into quasi-classical contribution 
and quantum corrections. The relative contribution of loops is characterized by the parameter 
ctgNc ex.p(ABFKLy), where y is the rapidity and Abfkl is the intercept of the BFKL pomeron 
(minus unity). So although the loop contributions are damped by factor 1/N% as compared with 
the tree contributions, which depend on a s N c , the large exponential factor makes their relative 
order grow with energy. It is of vital importance then to estimate the role of pomeron loops 
at present energies to have some idea of the validity of the quasi-classical methods which are 
currently applied to the analysis of the high-energy scattering in the framework of QCD, such as 
the Balitski-Kovchegov equation [2, 3], its generalization to nucleus- nucleus collisions [4] and 
the Color Glass Condensate approach (see the review [5] and references therein). 

In principle calculation of the loop contribution to the BFKL exchange in the lowest order 
seems to be straightforward. One has all the instruments for this goal: the BFKL propagator 
and the triple pomeron vertex. The main and only obstacle is the most complicated form of 
the latter, which converts a realistic calculation of the loop contributions into a formidable 
task. For these reasons there were quite a few attempts in this direction, all suffering from 
crude approximations and uncertainties. In [6, 7] the pomeron loop was calculated under the 
assumption that the triple-pomeron vertex is independent of conformal weights. Apart from 
that, simplified (asymptotic) expressions were used for the BFKL propagator into which the 
loop is inserted, in spite of the fact that with the loop growing much faster this is not allowed. 
The results in [6, 7] do not agree (probably due to differences in numerical coefficients). In any 
case in [7] it was found that the magnitude of the loop is so small that it gives no significant 
contribution up to extraordinary high energies (rapidities of the order of 40). Later in [8] a step 
forward was taken to take into account the dependence of the loop on inner conformal weights. 



This allowed to study the relative magnitude of the loop for the BFKL propagator with a fixed 
(dominant) conformal weight. This relative magnitude was found to be much larger: the loop 
begins to dominate already at rapidities of the order of 10—15. Still this result did not solve the 
real physical problem, the contribution of the loop to the scattering amplitude, which is obtained 
after integration over all conformal weights and requires knowledge of the triple pomeron vertex 
as a function of three conformal weights. 

Actually the formal expression for the triple pomeron vertex as a function of three conformal 
weights was derived by G.Korchemsky long ago [9]. Nevertheless its complexity has till now 
not been used for practical calculations. In this paper we try to exploit Korchemsky's formulas 
to actually calculate the contribution of the pomeron loop to the hadronic scattering amplitude 
in full rigour. This contribution turns out to be quite large and becomes comparable to the tree 
result at rapidities of the order of 10. 

We have to note that in literature one can find several claims to sum all loop contributions 
in the colour dipole approach or in the so-called reaction-diffusion formulation of the scattering 
mechanism [10, 11, 12, 13, 14]. However very crude approximations made from the start do not 
allow to consider these results even minimally reliable. 

The paper is organized as follows. We recall some basic elements of the conformal approach 
in the next section. In section 3 we use this results to calculate the single-loop contribution to 
the scattering amplitude with more or less arbitrary form of the impact factors. The last section 
presents our numerical results and discusses influence of the loop contribution on the scattering 
amplitude. 

2 Conformal invariant technique 

The conformal technique is based on Mobius (conformal) invariance of the BFKL approach. It 
uses expansion of the pomeron propagator and triple-pomeron vertex in orthogonal functions of 
the conformal basis. The Mobius invariance drastically simplifies this expansion. One can use 
these results and rewrite expression for any Feynman diagram in terms of conformal basis. In 
this case integrations over gluon coordinates transform into summations over conformal weights. 
The BFKL pomeron propagator satisfies the BFKL evolution equation: 

(— + HBFKijgv-yirurxf!^) = S(y - y)VfV 2 2 5^{n - fi)^(r 2 - f 2 ), (1) 

where y is the rapidity, 



- (lnp? + ]np2 + ]n r 2 l2 p\ + 
z \ p l p 2 



Hbfkl = ^(bip? + \np 2 2 + — \ur\ 2 p\ + — \ur\ 2 p\ - 4^(1) (2) 



is the BFKL Hamiltonian with a = a s N c /n and pi and n are transverse momenta and coordi- 
nates of the reggeized gluons, r^ = n — rj. 

The solution of this equation was found in terms of expansion in the conformal basis [15]: 

9y-y'(n,r 2 ;f 1 ,f 2 ) = ^2 E li (r 1 ,r 2 )Ep(f 1 ,f 2 )gh(y - y'), (3) 

formed by functions E^. In complex notations they are 

where \i = {n, u, ro} denotes a set of conformal weights n, v and two-dimensional center-of-mass 
coordinate vq. Integer n and real v enumerates functions of the basis. The conformal weight 



defined by n, v is 

1 4- n 

h= +iv; h = l-h*. (5) 

For simplicity we shall also use h to denote a set {n, v}. Restriction /j, > means using half of 
the whole set with u > 0, so that 
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§ = „5J« d V/ dV »" = T^T^7I- (6) 

Passing from rapidity y to the complex angular momentum j = 1 +w, we have the pomeron 
propagator 

/CO 
dye^ y g y (r 1 ,r 2 ;fi,f 2 ) (7) 

-co 

with the inverse transformation 

r-a+ico ^ 

gy(ri,r 2 ;fi,f 2 ) = -^e ujy g UJ (r 1 ,r2;f 1 ,f 2 ). (8) 

In fact — cj is the "energy" corresponding to the BFKL Hamiltonian. 
As a function of w, in the conformal basis, the propagator is 

9u,,h = t , (9) 

where w/i are the BFKL energy levels 

u h = 2a s (V(l) - Re^(h)) (10) 

and 

4tt 8 
//,= • (11) 
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The form of the triple-pomeron vertex V can be extracted from the interaction part of the 
Lagrangian of the effective non-local field theory [4], where pomerons are described by two 
non-local fields <3? and <I>^: 

2a 2 N c f d 2 r\d 2 r 2 d 2 r-i , , , + . . . , 

Li = — 2-^ / \ 2 2 - $(y, r 2 , r 3 $ y, r 3 , n)L 12 &{y, n,r 2 ) + h.c, 12 

n J r 2 2 r 2 3 r 2 3 

where operator L\ 2 is defined as 

L 12 = rf 2 V 2 V 2 (13) 

and acts on splitting pomeron or the one formed after fusion. Passing to the conformal repre- 
sentation one presents: 

T( ri ,r 2 \r 3 ,r 4 ;r 5 ,r 6 )= ^ r^^E^n,^)^ (r3,r4)£* 3 (r 5 ,r 6 ). (14) 

Mi>M2,M3>0 

The conformal invariance of the pomeron interaction dictates the structure of Lj^. as: 

r Mi|M2,M3 = i??2 12 ^23 3 ^3 13 >< (CC) X rfa, /*2, /*), (15) 

where parameters a^ are known function of conformal weights n, v, see [1], and (c.c.) refers 
to the complex conjugate of the preceding factors. 
In the lowest order 

Cl™ = ^2 12 #23 23 ^3 3 X (CC.) X Q(h U k2, h). (16) 



The conformal vertex fi was introduced and studied by Korchemsky [9] and corresponds to 
planar diagrams contribution, which gives the dominant part in the limit of large number of 
colours. As mentioned in the Introduction, in [9] explicit formulas for £1 were derived. In 
our numerical studies we use the results of Korchemsky in the following form (they are slightly 
different than in [9]): 

n(/n, h 2 , h 3 ) = 7r 3 [r 2 (/ Jl )r 2 (/ J2 )r(i - hi)v(i - h 2 )v(i - ha)]- 1 

3 

x ^2J a (hi,h 2 , h 3 )J a (hi,h 2 ,h 3 ), (17) 

a=l 

where functions J a and J a are defined via convolutions of hypergeometric functions: 

Ji(hiMM) = r(/*i + h 2 - h 3 )v{\ - / 11 )r(/ il )r(i - h 2 )v{h 2 ) 

x / dx (l-x)- h3 2 F 1 (hi,l-hi;l;x) 2 F 1 (h 2 ,l-h 2 ;l;x); 

Jo 

r(/M + h 2 - h 3 )T(i - / Il )r(/ Jl )r(i - h 2 )v{h 2 )v 2 {i - h 3 ) 



J 2 (hi,h 2 ,h 3 ) 



r(l + hi - h 3 )T{2 -fn-ka 
h 2 , 1 - h 2 , l-h 3 , 1 - h 3 
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and 



J 3 (hi,h 2 ,h 3 ) = J 2 (h 2 ,hi,h 3 ) (18) 

JiChiMM) = (-i) ni+n2 - 



T(-h! + h 2 + h 3 )v(i - h 2 )v(i - /n)r(/n; 



T(-hi + h 3 + 1) 
dx x~ hl (l - x) h2 ~ l 2 Fi{h 2 , -hi + h 2 + h 3 ; -hi + h 3 + 1; x) 2 Fi(l - h 3 , 1 - /h; l;x); 

1( - h - h - h , , 1 , W1 r(i ~ ^i)r(fe 3 )r(-fei + fr 2 + fe 3 )r(i - fe 2 ) 

J2(/tl,/l2,^3) = (-l) 2( r ^7 , u 

r 2 (-/ii + /i3 + i) 
x / dxx h3 - /ll (i-^) fe2_1 2i ? i(i-/ii,i-/ii;-/ii+/i3+i;^)2i ? i(/i2,-/ii+/i2+/i3;-/ii+/i3+i;^); 

JO 

•M^i, fa, fa) = ^(fa, fa, fa)- (19) 

In any pomeron Feynman diagram BFKL propagators and triple vertices are convoluted by 
integration over coordinates or momenta. One can perform this integration using expansions 
(3), (14) and completeness condition for the conformal basic functions: 

^p^E^n^E^n^) = V- (20) 
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As a result this convolution is substituted by the summation over conformal weights and inte- 
gration over intermediate center-of-mass coordinates. The latter can be done analytically using 
the known dependence (15). 

In the energetic representation all propagators become functions of w. The interaction vertex 
conserves the energy. Indeed, a product of two pomeron propagators, which form the simplest 
loop, contains the following integration over intermediate energy in w-space (note that integra- 
tion over gluon coordinates is suppressed): 

dye~ wy gi{y)g 2 (y) = / — £1 (u)i)g 2 (u - wi)- (21) 

JfT—ioc ^TVI 



It is obvious that number of integrals over energy corresponds to the number of loops of the 
Feynman diagram. 

With a help of this rules one can immediately write down the Schwinger-Dyson equation for 
a full conformal pomeron propagator G u h, which includes arbitrary number of loop insertions: 

Guj t h = 9u,h ~ 9u,h<>h*-'w,hS'u,h-i (22) 

where l\ comes from the triple-pomeron vertex operator L acting on the incoming and outgoing 
pomeron propagators. Function £ w 7, is the pomeron self-mass in the conformal basis: 

SW = 8 -^fi I ^ V r«? c w c w r w , 2 „. (23) 

M1M2 

Using (15) one can perform integrations over intermediate center-of-mass coordinate explicitly. 
As a result we are left with 



_ 8aJN? f dui ^ (o) . . 

^", h - n 2 J ~^~ 2^ l h\hiM u ' hl U "M L h u h 2 \h, W 



8aiN 2 [duj. 

i 

hih.2 

where 



T. = i t r *>(* + $)■ (») 



7T — jn 
n=— 00 " 



The formal solution of equation (22) is 

Gu > h = TTn — TJiy — • ( 26 ) 

In the lowest order one should change G w ,h — >• 5W and r^^ — >• T^ fe , h in (24). 

3 Scattering amplitude 

The amplitude for the scattering of two hadrons at rapidity y can be presented as 

rcr+ioo A,. 

A(s,t) = is ^S^Uq 2 ), (27) 

where function fu(q 2 ) is a t-channel partial wave for fixed transferred momentum. 

In the approximation of a single bare pomeron exchange it is a convolution of the bare 
pomeron propagator g u with two impact factors describing interaction of the pomeron with 
external hadrons. Passing to the coordinate space we have 

(2ir) 2 fUq 2 )S {2) (q-q')= I d 2 nd 2 r 2 d 2 r 3 d 2 r 4 ^i{n,r 2 ,q)g w {n,r 2 \n,r A )^* 2 {r^r i ,q t ). (28) 



The impact factors ®i(ri,r 2 ,q) in the coordinate space in (28) are related to those in the 
momentum space by 

Hn,r 2 , q) = e* R *{r, q), Mr, q) = J ^ *(*, q ~ k)e^~i>, (29) 

where R = (n +r 2 )/2 and r = r\ 2 . Note that for a colorless object the impact factor &(k, q — k) 
in the momentum space have to vanish at k = 0. Correspondingly in the coordinate space 
integration of <J>(r, q) over r has to give zero. This can always be achieved by changing 

$(r, q) -)• $(r, q) - 5 2 (r) f d 2 r$(r, q). (30) 



However the term with the (5-function does not contribute, since the pomeron propagator in (28) 
vanishes at r\i = or r^ = 0. So we are not to bother about this property of <&(r, q) for a 
colorless target provided it is integrable over r. 

Putting (29) into (28) and separating the dependence on center-of-mass coordinates we get 



Uq 2 ) = j d 2 rd 2 r' <h(r, q)&{r, r')<S>*(r', q), (31) 

where the pomeron propagator in the mixed representation is 

(27r) 2 gUr,r')5^(q - q') = f d 2 Rd 2 R' e^ R e^' R ' g^n, r 2 |r 3 , r 4 ). (32) 



Jui,h 



Here g^(r,r') describes propagation of the pomeron with momentum transfer q and can be 
interpreted as the amplitude for the scattering of two dipoles with sizes r and r' . 

The expression g^(r,r') in the conformal basis can be easily found if one uses the corre- 
sponding conformal expansion for the bare propagator. Indeed the pomeron Green function in 
the conformal basis is 

2 r roo , n 2 \ 

gu(ri,r 2 \r 3 ,r 4 ) = —^ ^ / d 2 r / dv\v 2 + — )gu,hE^(r w ,r 2 o)E^(r 30 ,r m ), (33) 

so that 

x j d 2 R e^ R E,(R+ r -,R- r -) j d 2 R' e~^ R ' E;(R' + r -,R' - r -). (34) 

For our purposes it will be sufficient co consider scattering at zero transferred momentum, 
i.e. forward scattering amplitude /^(0). For this we have to know behavior of the last two 
integrals in (34) at q 2 — >■ 0. The Fourier transformation of conformal functions (4) in the region 
of small momentum transfer was investigated in [15]. One can use this results to pass to the 
limit q = 0: 

tiM = - 2 \rr'\j: I du\-\ (— ) ' g^ h . (35) 

We shall be interested only in the leading contribution, which comes from n = 0. With n = 

9°Jr,r') = - / dv |r| 1+2i >'| 1_2il W (36) 

7T Jo 

Therefore the forward scattering amplitude is 

/«(°) = ^2 J^ dv 9^u J d 2 r$(r,0)\r\ 1+2w J d 2 r'<5>*(r' Mr 1 ^ 2 ™ . (37) 

It is natural to choose the impact factor in a Gaussian form 

$(r,0) = — e~ br \ (38) 

IT 

where b is the inverse of the hadron radius squared b = l/R% and dimensionless A is the effective 
coupling of the hadron to pomeron. The value of A is obviously irrelevant for the study of the 
relative contribution of loops. With choice (38) integration over r can be trivially done: 

[ d 2 r<S>(r,0)\r\ 1+2iu = — Tr^ 3 / 2 ^^ + it/). (39) 



Putting this result into (37) we find 

1 f'°° A 2 / 1\ 7T 

fU0) = — dug UiV x— (v 2 + -)—— — (40) 

ir z Jo b V 4/ cosn(7r^j 

where <7 W ^ is a conformal propagator (9) with n = 0. 

To take into account the loop contribution one has to substitute the bare propagator g uv 
by the full Green function. With a single loop insertion 

r - - —— Sa ^ (£\\ 

!/&>,«/ + Iq u ^oj,u 'Oi/ ^ - UJ V (W - W„)^ 

where now 

w„ = 2^(^(1) -Re ^(1/2 + «/)) (42) 

and 

1 1 1 



W 16(^ + 1/4)2- (43) 

The first term in (41) comes from exchange of the bare pomeron. It is not difficult to show 
that in terms of rapidity y it gives 

A^(0) = -^— r dv -. l . -^— e^. (44) 

y K ' 16vr 2 b Jo ( v 2 + l\ cosh(7rz/) V ' 

The second term in (41) corresponds to the lowest order loop contribution. As a function of uj 

f<?)(0) = l —^- r dv idJ + \) , S "'\ 9 *, , (45) 

Ju v J 16vr 2 b Jo V V (u - u„) 2 cosb(irv) v ' 

with the explicit form of X Wj „ given in [8]: 

a t N c f°° j,. j,. _^i A_ n 2 (l/2 + «i/, 1/2 + ii/i, 1/2 + iiAj 



X^ = „ s in c / dv\dv2- -^- -2 77; r 77; r (46) 

8vr 10 Jo ( v 2 + iV (2 + 1 V w-w(0,i^)-w(0,iA2) 



and the conformal vertex given by (18) and (19). 

For numerical calculations it is convenient to pass to rapidity. Performing integration over 

co we find 

1 A 2 r°° , / o 1\ tt 



^•)=-is?Tr*" 8 ^ + i) 



COsh(7T^) 



\ dv x dv 2 ? n 2 (l/2 + ii/, 1/2 + iv x , 1/2 + «/ 2 ) 

•/0 f,,2 1 l\ /\,2 1 l\ 



87r10 7 ° ^ + i) \»2 + I 

" 2+ 2 )• ( 47 ) 

w„ - o;^ - w^ ( w „ - Wvi -uv 2 ) (wj, — u vi - uj V2 ) ' 

The total forward scattering amplitude with the lowest order loop correction is a sum of (44) 
and (47): 

^(0)=4 1 )(0)+4 2 )(0). (48) 

Numerical calculations of (44) and (47) will be presented in the next section. 



4 Numerical studies 

We have set up a program which calculates the bare pomeron exchange amplitude (47) and single- 
loop contribution (44) . By far the most difficult part is the computation of the triple-pomeron 
vertex Q. It requires complicated numerical procedures and is extremely time-consuming. The 
vertex is a complex function which depends on three conformal variables v. We have restricted 
this variables to lie in the interval < v < 3.0 and introduced a grid dividing this interval into N 
points. We found the vertex on the grid, which requires calculation at iV 3 points. This strongly 
limits numbers N admissible for given calculation resources. In our case we used N = 20 
compatible with reasonable computation time. The value of the vertex in between the grid 
points was found by interpolation. 

The integrals in (44) and (47) were calculated by the Newton-Cotes integration formulas. 
The limits in v were taken as before < v < 3.0. The number of sample points was chosen to 
provide relative error 10 -3 . The same integration strategy was used for the vertex integrals in 
(18), (19). 

We have performed calculations for the standard value of the QCD coupling constant a s = 0.2 
and N c = 3. In Fig. 1 the dash-dotted line presents the bare pomeron exchange contribution 
(44). The behavior of the amplitude is determined by the initial pole of the conformal BFKL 
propagator (9). As expected, the curve grows with rapidity roughly as e y , where A ~ 0.48. 
The behavior of the single-loop contribution (47) is shown in Fig. 1 by the solid line. It roughly 
grows twice faster, as ~ e y , again as expected. Note that the results in Fig. 1 are normalized 
by factor 16vr 2 6/A 2 . 

For small rapidities the loop term is suppressed by the smallness of the QCD coupling 
constant. However, its faster growth with rapidity compensates this very early. This is illustrated 
in Fig. 2, where the ratio R(y) = \A y (0)/Ay (0)| is shown (recall that this quantity does not 
depend on the hadron radius nor effective coupling of the hadron to pomeron A). One clearly 
observes the exponential growth of R: R(y) ~ exp(Ay). 

As a result we conclude that the loop contribution becomes visible already at rapidities 3 — 8 
and starts to dominate at y ~ 8 — 10. 

5 Conclusions 

We have studied a single loop contribution to the scattering amplitude of two colliding hadrons. 
We have found expression for the amplitude in a framework of conformal invariant technique with 
more or less general form of the impact factors. The triple-pomeron vertex with full dependence 
on the intermediate conformal weights was calculated and used. 

Numerical analysis shows that smallness of the QCD coupling constant is compensated by 
rapid growth of the single-loop amplitude with rapidity. We found that the loop contribution 
manifests itself at relatively small rapidities y ~ 3 — 8 and dominates the bare pomeron exchange 
amplitude already at y ~ 8 — 10. In general this conclusions are in accordance with [8]. 

Thus loops have to be taken into account already at present energies. Higher order loops 
calculation is required for larger values of rapidity with summation of all loops in the limit 
y — > oo. This nontrivial problem is left for future investigations. 
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Figure 1: The bare pomeron exchange (44) (dash-dotted curve) and singledoop (47) (solid curve) 
contributions to the forward scattering amplitude as functions of rapidity. 




Figure 2: The ratio \A y ' (0)/Ay (0)\ as a function of rapidity. 



References 

[1] M.A.Braun, Eur. Phys. J C 48 (2006) 511. 

[2] I.Balitsky, Nucl. Phys. B 463 (1996) 99. 

[3] Yu.Kovchegov, Phys. Rev. D 60 (1999) 034008; D 61 (2000) 074018. 

[4] M.A.Braun, Phys. Lett. B 483 (2000) 115; Phys. Lett. B 632 (2006) 297. 

[5] F.Gelis, E.Iancu, J.Jamal-Marian an R.Venugopalan, Ann. Rev. Nucl. Part. Sci 60 (2010) 
483. 

[6] R.Peschanski, Phys. Lett. B 409 (1997) 491 

[7] J.Bartels, M.Ryskin and G.P. Vacca, Eur. Phys. J. C 27(2003) 101. 

[8] M.A. Braun, Eur. Phys. J. C 63, 287 (2009); 

[9] G.Korchemsky, Nucl. Phys. B 550 (1999) 397. 

[10] J.-P. Blaizot, E.Iancu, K.Itakura and D.N.Triantafyllopulos, Phys. Lett. B 615 (2005) 221. 

[11] E. Levin and M.Lublinsky, Nucl. Phys. A 763 (2005) 172 

[12] E.Iancu and D.Triantafyllopulos, Phys. Lett. 610 (2005) 172. 

[13] S.Bondarenko, Nucl. Phys. A 792 (2007) 264 

[14] E. Levin, J.Miller and A.Prygarin, Nucl. Phys. A 806 (2008) 245. 

[15] L.N. Lipatov, Zh. Eksp. Teor. Fiz. 90 1536 (1986); in " Perturbative QCD", ed. A.H.Mueller, 
World. Sci. Singapore (1989). 



10 



